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Abstract 

Charge carriers in bilayer graphene occupy two parabolic continua of electron-like and hole-like 
states which differ by the alignment between carrier pseudospin and its momentum, the property 
known as chirality. Due to chirality conservation, a strong confining potential can host unusual 
bound states: electron levels cloaked into the hole continuum. The energy levels and the wave 
functions of the cloaked states can be obtained by solving the Schrodinger equation for a massive 
non-chiral particle in the p-wave channel in two dimensions. Eventually, cloaked states slowly decay 
into the continuum, via trigonal warping effects. We discuss the key properties of cloaked states 
in circularly symmetric potentials, and show that cloaking should be observable in quantum corral 
geometries via scanning tunneling probe measurements. 
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There are two related fundamental facts that make the electronic properties of graphene 
interesting: the coexistence of electron and hole states in the vicinity of conic points in the 
spectrum, and chirality of charge carriers. Charge carriers in graphene can be characterized 
by their pseudospin, and different alignments between pseudospin and momentum result 
in electron-like and hole-like states. There is an important difference between the cases of 
monolayer and bilayer graphene (BLG) due to their different dispersion laws and different 
structure of chiral states. In single-layer graphene, the linear dispersion results in the density 
of states vanishing at the Dirac point [T]. In bilayer graphene, the dispersion is parabolic, and 
the density of states is flat [2] , so that there is a large number of both electron- and hole-like 
states near the point separating the two subbands. One consequence of this is a potential 
instability of bilayer graphene with respect to formation of broken-symmetry phases j3]. The 
chirality also affects transport properties. In monolayer graphene, any potential barrier is 
perfectly transparent at normal incidence [1]. In bilayer graphene, chirality is responsible 
for electron cloaking [3] described below. 

The effects of chirality are manifested in the properties of bound and quasi-bound states. 
In monolayer graphene, a strong conhning potential, such as the supercritical Coulomb 
potential[6], can give rise to a quasi-localized resonance which was recently observed in 
scanning tunneling microscopy (STM) experiments [7]. The resonant state is formed by 
electron-like states submerged into the hole band and coupled to it via Klein tunneling. The 
number of states available for decay is suppressed by linear density of states, so that the 
width of the resonant peak is smaller than its energy, and the resonance can be resolved. 

Do localized resonances exist in BLG? At hrst glance, the large flat density of states in 
BLG would result in the fast decay of such a state. However, the recent analysis of of barrier 
transmission problem in bilayer graphene |5] suggests that the coupling to the continuum can 
be suppressed as a consequence of chirality conservation. For a one-dimensional barrier U{x), 
the quasimomentum perpendicular to the barrier, py, is conserved. For normal incidence, 
Py = 0, electron and hole continua decouple, and the problem is equivalent to two one¬ 
dimensional Schrodinger equations for two massive, non-chiral particles, which differ only 
by the sign of the potential energy. An attractive potential U{x) < 0 hosts at least one 
bound state at negative energies, for any potential strength. The state is submerged into 
the continuum of holes, repelled by the potential —U{x). In |5] this phenomenon was referred 
to as cloaking, by analogy with cloaking in optics [8l|9]. Indeed, one cannot probe, e.g., the 
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occupancy of the level using decoupled hole states coexisting in energy. Away from normal 
incidence, however, the coupling between electron and hole states is proportional to so 
that the localized level acquires hnite width oc p^, and thus is turned into a resonant barrier 
state which can contribute to the transmission [10]. Thus, a cloaked state appears as a 
singular point at Py = 0. Another manifestation of decoupling between electrons and holes 
in quantum dots in graphene was mentioned in the context of quantum dots in BLG HU. 

In this Letter, we show that cloaked states in BLG can be hosted by a local, circularly 
symmetric potential. The cloaking occurs in the s-wave channel, which serves as an analogue 
of the normal incidence channel in the barrier transmission problem. Thus, the cloaking 
relies on a discrete quantum number, which results in perfect isolation of the cloaked state 
from the continuum, and gives a truly discrete level. Unlike the case of one dimension, 
however, we also hnd that the bound states only occur if the potential exceeds a certain 
critical strength. We also consider how the cloaked states decay due to trigonal warping, and 
show that their resonant width is less than ImeV. We also discuss potential experimental 
signatures of cloaked states. 

To analyze the dynamics of charge carriers in BLG, we introduce the following single¬ 
particle Hamiltonian [2] 


H = 


1 

2m* 


0 iPx - iPyf 
(Px + iPyY 0 


+ U{x,y) , 


( 1 ) 


where Px^y = —ihdx^y are the standard momentum operators, m* O.OSbrrie is the effective 
mass of charge carriers, and U{x,y) is the conhning potential. The matrix Hamiltonian 
operates on two-component wave functions which represent the two pseudospin states. For 
free electrons {U = 0), the eigenstates of this Hamiltonian give two parabolic dispersion 
bands, Ck = ±h^k^/(2m*). The chirality of charge carriers is reflected in the non-trivial 
phase of the momentum combination (p^, -t- ipyY- In particular, this leads to the relative 
momentum-dependent phase between the two components of a plane wave: T 2 = 
where ^k is the azimuthal angle of the wavevector k. 

We analyze the case of a circularly symmetric potential, introducing the polar coordinates: 
X = rcos0, y = rsin0, so that U = U{r). The quantum states can be characterized by 
their angular momentum quantum number M = ..., —2, —1, 0,1, 2,..., which takes integer 
values. We shall represent the wave function as a superposition of two components of 
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opposite chirality described by the two amplitudes, um{j) and vm{j)'- 


d'Ar(r,0) = e 


_ iM(l> 


( 2 ) 


e ^‘>'[uM{r)+VM{r)] 
e"’^[uM{r) -VM{,r)] 

Note that the magnetic number for the two pseudospin components differ by two; this 
is characteristic of bilayer graphene. Let us substitute the ansatz (|^ into Schrodinger 
equation H'^m = m with the Hamiltonian ([^. In polar coordinates, the differential 
operators dx ± idy take the form [dr ± • The chiral components UM{r) and vuir) 

can then be separated by calculating the sum and difference of the two equations: 


[e -U]um = 
-[e - U]vm = 




2m* 

JL 

2m* 


T^mUm — 
'^mVm — 


2Mv'n 


M 


r 

2Mu'f 


(3) 


M 


r 


- „ u — 1 . du 

where Vmu = u H-1- -—u ,u = — . 

r dr 

The terms proportional to 2M/r describe the coupling between umI'^) and and makes 

this system equivalent to a biharmonic equation. Interestingly, the equations decouple in 
the s-wave channel, M = 0: 


^“0 = + U{r)uo , ( 4 ) 

2m 

“^^0 = - U{r)vQ . 

2m 

These equations resemble the ones considered in [5], in the limit of normal incidence. The 
two chirality components are described by the two decoupled amplitudes Uo{r) and Vo{r); the 
conhning potential U (r) is attractive for one component and is repulsive for the other. The 
differential operator Vq on the right-hand side of Eqs.(|^ is identical to the radial part of the 
Laplacian operator in two dimensions in the angular momentum channel with M = ±1, i.e., 
in the p-channel. Therefore, the solutions in the s-wave channel for bilayer graphene could 
be constructed from p-wave solutions of the Schrodinger equation for a massive non-chiral 
particle in two dimensions. (Alternatively, one can map the wave function onto the s-state 
in four dimensions, by the change of variable Mo(^) = '^WQ{r), which brings the kinetic energy 
to the form Wq +3wQ/r.) For positive chirality states described by uo{r) the potential in the 
Schrodinger equation is U{r), and the energy is e; for states in the opposite chirality channel 
these quantities change sign. Thus, the two components, uo{r) and vo(r) describe the two 
s-wave subbands of electron-like and hole-like states, respectively. 
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Before employing this mapping, let us recall the properties of bound states for a massive 
non-chiral particle. In dimensions less than or equal to two, a bound state in the s-channel 
exists for an arbitrary weak potential; however, the binding energy decreases with increasing 
dimension na. Bound states in higher momentum channels exist only if the potential is 
sufficiently strong, so that it can overcome quantum zero-point motion energy and the 
centrifugal barrier. The massive particle in more than two dimensions can be conhned by 
a negative potential U (r) if its strength exceeds a certain critical value, of the order of 
the zero-point energy in the potential. For a three-dimensional square well [13] of radius a 
(f/(r) = —Uq for r < a), this value is Uq^ct = ^rn*a? • same analysis can be carried over 
for a bound state in the p-channel in two dimensions: the bound state becomes detached 
from the continuum when the wave function inside the well reaches its maximum at the well 
boundary. On the other hand, the wave function in the p-channel regular at r = 0 must 
be proportional to r. For a rectangular well, the wave function inside the well is Ji{qr)^ 
where q = ^2m*f/o/h is the wavenumber inside the well, and Ji{x) is the hrst order Bessel 
function. This function reaches its hrst maximum at a; = 1.84, which gives the critical well 
depth t/o,cr = 1.70h^/(m*a^). The estimate ffo,cr ~ h?/{m*a?) is valid irrespectively of the 
details of the well prohle. 

We can now state the central hnding of this paper: discrete bound states exist in the 
electron spectrum of bilayer graphene when the potential U{r) exceeds the critical strength: 
if it hosts a bound state in the p-channel for the massive non-chiral particle in two dimensions, 
it can also host a cloaked bound state in BLG. The critical potential strength is of the order 
of the kinetic energy of quantum zero-point motion, h?/{m*a^), where a is the radius of the 
potential. Note, however, that the bound state in bilayer graphene occurs irrespective of the 
sign of the potential energy, due to the particle-hole symmetry. For an attractive potential, 
the cloaked state is split off the electron continuum and submerged into the hole continuum, 
while a repulsive potential gives rise to a hole-like state. 

Let us now discuss the physical situations in which the cloaked states can occur. It 
is very natural to check if cloaked states can be hosted by charged impurities. One can 
consider the cloaking by the Coulomb potential U{r) = —Ze^/ xr, induced by the charge Ze, 
where x is the effective dielectric constant of graphene and its environment. For such a 
potential, inhnitely many Rydberg-like bound states exist at arbitrarily small charge Z in 
all momentum channels. In principle, a Coulomb potential acting on electrons in bilayer 
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FIG. 1. The position of the cloaked level as a function of the total charge for the two distribu¬ 
tions: (i) point-like charge (black), (ii) charged ring (red, blue). One can see that the cloaking 
is facilitated by deconcentrating the charge. The screening is slightly affected by the background 
charge density n^: for the lower charge density (the blue curve), the screening is enhanced at 
short distances owing to the wavevector dependence of the polarization operator m- This re¬ 
sults in smaller energy of the cloaked state. The dashed lines indicate the position of the Fermi 
level. Inset: the potential energy profile U{r) in the ring configuration for the background charge 
density n* = 3 • 10^^cm“^. 


graphene would also host a hydrogen-like spectrum of cloaked states na, with energies = 
—ErZ^I{ n+3/2)^, with Rydberg energy Er = \m*e^/ (xh)^. For bilayer graphene on a Si02 
substrate (x = 2.5), the energy of the lowest cloaked state would be Cq = —35meV, due to 
small effective mass. However, unlike the case of donor states in semiconductors, an external 
Coulomb potential is effectively screened in bilayer graphene due to the non-vanishing density 
of states and = 4 fermion flavours. The corresponding Thomas-Fermi screening radius, 
h?>c/{Nm*e^) ~ 0.37nm is therefore shorter than the Bohr radius h^x/(m*e^) ~ 1.5nm. 
The screening by charge carriers in BLG is mostly similar to Thomas-Fermi screening [T5] . 
except at distances shorter than the Fermi wavelength, where the response is larger by the 
factor of log 4 1.4. This response function can be employed to hnd the potential induced 

by a point charge, and then solve Eq. (|^ for uo{r) numerically for different values of the 
Coulomb charge Z. The resulting dependence of the cloaked level position on Z is shown as 
the black curve in Fig. [Tj The cloaked state exists if Z exceeds the critical value Zc ~ 12. 
Therefore, it is very unlikely that cloaked states arise as donor or acceptor states due to 
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FIG. 2. The Local density of electronic states (LDOS) in quantum corral formed by a ring 
of charged impurities (3nm in radius, total charge Z = 12, background electron density = 
3 • lO^^crn”^) was obtained by solving Eq.([^ numerically for all channels with |M| < 10. The 
trigonal warping contributions were included for the M = 0,3, —3 channels. The LDOS is shown 
in the units of its bulk value, vq = m*/{2TTh‘^). Near the resonant energy, the LDOS exhibits a 
sharp peak at e = —86meV that represents the probability distribution in the cloaked state. The 
peak dominates over the surrounding continuum due to its small width, 0.64meV. Inset: the LDOS 
as a function of energy and the distance away from the center of the corral. The resonant peak is 
seen as a sharp streak superimposed onto the smooth interference fringes formed by incident and 
repelled waves in the hole continuum. 

charged impurities. The culprit is the repulsion by a centrifugal barrier in the p-state, which 
prevents the particle from reaching the attractive core of the potential. 

The estimate for the critical potential strength, Ucr ~ h?/{m*a‘^), shows that cloaked 
states can arise in smooth potentials with a large characteristic distance a. For example, 
if a = lOnm, the critical strength is estimated to be of the order of SOmeV. One can also 
notice that distributing the positive charge rather than concentrating it at a single point is 
a way to get around the centrifugal barrier. This suggests a better way to engineer cloaked 
states. One can create a large cluster of charge impurities by manipulating them with the tip 
of a scanning microscope. Arranging the impurities into a charged ring (such as a quantum 
corral [IE]), one can lower the required threshold. For example, for impurities placed around 
the ring of 3nm radius, the simulations give a lower value of the critical charge, Zc = 7, 
weakly dependent on the background electron density, as shown by the red and blue curves 
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in Fig. It should be also possible to engineer cloaked states via gating a BLG sample 
locally, e.g, with the tip of a scanning tunneling microscope [HI- 

To detect the cloaked resonance, one can measure the local density of states (LDOS) via 
scanning tunneling microscopy measurements. The LDOS can be dehned in terms of the 
wave functions {uM,e{r),VM,e(j')) at a given energy e as z^(e, r) oc sumM[u %[+ 

To show how the bound state could reveal itself in STM measurements, we have calculated 
the LDOS numerically, including the effects of RPA screening. The decay due to trigonal 
warping (see below) was included as Lorentzian broadening of the level. Figure]^ shows the 
LDOS in the vicinity of a charged ring with total charge Z = 12. The resonant contribution 
can be seen as a sharp £n at Cc ~ —86meV sticking out of the hole continuum. 

We now turn to the issue of stability of cloaked states. Deviations of the kinetic energy 
from the chiral form ([^ results in the decay of the bound states. The most important 
next-order contribution to the kinetic energy is the trigonal warping term [2]: 


H-u] ^3 


0 Px + ipy 
Px - ipy 0 


(5) 


where the parameter ua ~ 6 ■ lO^m/s. The warping breaks chirality and angular momentum 
conservation. Applying to the wave function (|^, we hnd 





Mo(r) 


—ihvs 



r - 



( 6 ) 


Hence, the cloaked M = 0 state could decay into M = ±3 continuum states. This decay, 
however, is partially suppressed by two effects. Firstly, the wave function with large angular 
momentum is suppressed near the origin r = 0 due to the centrifugal barrier. Second, the 
continuum state, unlike the cloaked state, is repelled by the potential —U{r). Near the 
critical threshold, when the resonant state is shallow, the dominant effect is due to the 
centrifugal potential. A simple estimate of the resulting resonant width can be found with 
the help of Fermi’s golden rule. Let us assume, for simplicity, that we consider the state near 
the cloaking threshold. This state can be described by a shallow wave function, independent 
of the prohle of the potential. For the p-channel, the decaying shallow-state wave function 
is proportional to the hrst-order Macdonald function Ki{Kr). The parameter k is related 

fc2 2 

to the energy of the cloaked state: Cc = — 2 A* outgoing wave that coexists at the same 
energy has the wavenumber k, and the corresponding wave function in the M = 3 channel 





can be written in terms of Bessel fnnctions The properly normalized initial and 

final states are 


KKi{Kr) / e \ 


Tout = 



I T2(Kr)e^*'^ 
y J4(Kr)e'^*'^ 


(7) 


where i? is a large normalization radins for continnnm states. Applying recnrsive relations 
between Bessel fnnctions and known integrals [H], one can hnd the transition matrix ele¬ 
ment: 

1/ = <^d^out|i7u.|Tin) = -thvs^ 16i?log A • 

The density of hnal states in a given angnlar momentnm channel is i'm =±3 = (-R/vr) {dn/dE) = 
{Rm*)/{ ttIi?k). The golden rnle formnla then gives the decay width as 


O^r TT 

T = 2x-\V\^VM=3 = ^^rrfvl (9) 

(we have inclnded both momentnm channels here). Thns, we see that the broadening 
of the cloaked state is only weakly dependent of its energy, and is given by the qnan- 
tity (7r/4)m*n| 0.6meV. Therefore, the warping becomes insignihcant when the energy of 

the cloaked state exceeds this valne. By dimensional argnments, this estimate shonld hold, 
by the order of magnitnde, irrespective of the details of the well U{r). 

One shonld also bear in mind that transverse electric helds can open a gap in bilayer 
graphene [I9]. The gap A violates chirality conservation, and gives the resonance a hnite 
width, T ~ A^/cc; the resonance is qnenched if the gap is larger than its energy ec- Thns, the 
resonance can be controlled by a transverse held, and this can be potentially employed to 
identify the signatnre of cloaked resonances in the experimental data. The cloaked state can 
be also destroyed by magnetic held. Indeed, the ehect of the magnetic held B amonnts to 
shifting the angnlar momentnm. The resonance is destroyed when the magnetic hnx throngh 
the cloaked orbit becomes of the order of one hnx qnantnm, or when the cyclotron energy 
is comparable with the energy of the cloaked state. This analysis also shows that it may 
be desirable to avoid non-nniform strain in the cloaked region to minimize the detrimental 
ehects of psendomagnetic helds [20] . 
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An interesting open question is the role of symmetry in this problem. Cloaked states 
have been found here for strong potentials with rotational symmetry; earlier, cloaked states 
were identihed for an arbitrarily weak potentials with translational symmetry [S]. The two 
solutions share a common feature: the two components of the wave function of the cloaked 
states are complex conjugate to each other: T = ('0(a^, 2/), y))- Indeed, if the Hamil¬ 

tonian admits solutions of this type, one can show that the two chirality components are 
decoupled. Whether this can occur in other geometries, in particular, non-symmetric ones, 
perhaps, with a different potential strength threshold, is not known to the author. 

Thus, we have shown that chirality of quasiparticles in bilayer graphene results in un¬ 
conventional bound states that do not hybridize with the surrounding hole continuum. 
Such states occur if the conhning potential exceeds the quantum zero-point motion en¬ 
ergy: \U\ > The cloaking can be implemented with both attractive and repulsive 

potentials. Trigonal warping converts these states into resonances of hnite width, transverse 
electric and magnetic helds can be employed to quench these states. We predict that such 
states could be engineered in large clusters of charged impurities, and probed by scanning 
tunneling microscopy experiments. The author would like to thank L.S.Levitov, V.I.Falko, 
and P.Brouwer for useful discussions. This research was supported by EPSRC/HEFCE No. 
EP/G036101. 
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